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VI}ration-induced  phase  noise  can  change  the  probability  of  detection  of  moving-target  indicator  (MTI)  radar  from 
near-100  percent  to  zero.  Oscillators  that  are  capable  of  meeting  the  requirements  of  MTI  radar  systems  in  a  quiet 
environment  are  readily  available.  In  the  vforati^  environments  of  airborne  platforms,  however,  the  phase  noise 
of  oscillators  degrades  very  signiTicantly.  For  example,  a  10-MHz  crystal  oscillator  may  have  a  phase  noise  of  -140 
dBc/Hz,  100  Hz  from  the  carrier.  Assuming  a  typinl  one  part  per  billion  per  g  acceleration  sensitivity,  this  phase 
noise  degrades  to  -93  dBc/Hz  under  a  0.1  sq  g^  random  vbration  at  100  Hz  away  from  the  carrier.  If  there  is  a 
1g  sinusoidal  vbration  at  100  Hz,  there  wiB  be  a  pair  of  -86  dBc  spectral  lines  superimposed  on  the  signal.  Of 
course,  upon  multiplication  to  10  GHz,  the  phase  noise  increases  by  at  least  another  60  dB. 

This  report  reviews  the  causes  and  effects  of  acceleration  sensitivity  of  bulvwave  quartz  crystal  resonators,  and  the 
methods  that  reduce  or  compensate  for  that  sensitivity.  Most  of  what  is  discussed  is  equally  relevant  to  surtace 
acoustic  wave  (SAW),  shallow  bulk  acoustic  wave  (SBAW),  dielectric  resonator-oscillators(DRO),  and  other  types 
of  oscWators. 
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INTRODUCTION 

Vibration-induiied  phase  noise  can  change  the 
probability  of  detection  of  moving*target  indica* 
tor  (MTI)  radar  'rom  near-100  percent  to  zero. 
Oscillators  that  are  capable  of  meeting  the 
requirements  of  MTI  radar  systems  in  a  quiet 
environment  are  readily  available.  In  the 
vibrating  environn^ents  of  airborne  platforms, 
however,  the  phsu'e  noise  of  oscillators  de¬ 
grades  very  significantly.  For  example,  a  10 
MHz  crystal  oscillator  may  have  a  phase  noise 
of  *140  dBc/Hz,  100  Hz  from  the  carrier,  where 
dBc  refers  to  dB  relative  to  the  carrier.  As¬ 
suming  a  typical  1x10**  per  g  acceleration 
sensitivity,  this  phase  noise  degrades  to  -93 
dBc/Hz  under  a  0.1  gVhz  random  vibration  100 
Hz  away  from  the  carrier.  If  there  is  a  1g 
sinusoidal  vibration  at  100  Hz,  there  will  be  a 
pair  of  -86  dBc  spectral  lines  superimposed  on 
the  signal.  Of  course,  upon  multiplication  to  10 
GHz,  the  phase  noise  increases  by  at  least 
another  60  dB. 

There  has  been  an  awareness  of  acceleration 
effects  in  frequency  sources  at  least  since  the 
advent  of  missile  and  satellite  applications,''* 
Doppler  radars,^**  and  other  systems  requiring 
extremely  low  noise.*  '®  There  has  not,  howev¬ 
er,  been  a  general  appreciation  of  the  magni¬ 
tude  of  the  problem.  As  a  case  in  point,  there 
is  little  or  no  mention  of  acceleration  sensitivity 
in  the  general  texts  on  quartz  crystal  oscillators 
and  resonators  or  of  vibration-ind'jced  phase 
noise  in  radar  textbooks. 

High-stability  frequency  sources,  including 
atomic  standards,  contain  quartz  crystal  reso¬ 
nators.  One  result  of  the  evolution  of  electron¬ 
ics,  i.e.,  the  transition  from  tubes  to  transistors, 
and  from  point-to-point  wiring  to  printed  circuits, 
is  the  establishment  of  the  qua.'tz  cryst^'l  reso¬ 
nator  as  the  most  acceleration-sensitive  com¬ 
ponent  in  frequency  sources.  This  ,>aper 
reviews  the  causes  and  effects  of  acceleration 
sensitivity  of  bulkwave  quartz  crystal  re&nna- 
tors,  and  the  methods  that  reduce  or  compen¬ 


sate  for  that  sensitivity.  Most  of  what  is  dis¬ 
cussed  is  equally  relevant  to  surface  acoustic 
wave  (SAW),  shallow  bulk  acoustic  wave 
(SBAW),  dielectric  resonator-oscillators  (DRO), 
atomic  frequency  standards,  and  other  types  of 
oscillators. 


THE  EFFECT  OF  ACCELERATION  ON  A 
CRYSTAL  RESONATOR 

A  quartz  crystal  resonator  subject  to  a  steady 
acceleration  has  a  slightly  different  series  reso¬ 
nant  frequency  than  the  same  resonator  exper¬ 
iencing  zero  acceleration.'  Furthermore,  it  has 
been  observed  that  the  magnitude  of  the  fre¬ 
quency  shift  is  proportional  to  the  magnitude  of 
ti.e  acceleration,  and  is  also  dependent  upon 
the  direction  of  the  acceleration  relative  to  a 
coordinate  system  fixed  to  the  resonator."  It 
has  been  shown,  empirically,  that  the  accelera¬ 
tion  sensitivity  of  a  quartz  crystal  resonator  is  a 
vector  quantity.'*  Therefore,  the  frequency 
during  acceleration  can  be  written  as  a  function 
of  the  scalar  product  of  two  vectors 

/(a)-4(1.f -a)  (’) 

where  f(d)  is  the  resonant  frequency  of  the 
resonator  experiencing  acceleration  S,  fo  is  the 
frequency  with  no  acceleration  (referred  to  as 
the  carrier  frequency),  and  T is  the  acceleration 
sensitivity  vector.  It  can  be  seen  from  Eq.  1 
that  the  frequency  of  an  accelerating  resonator 
is  a  maximum  when  the  acceleration  is  parallel 
to  the  acceleration  sensitivity  vector;  it  is  a 
minimum  when  the  acceleration  is  antiparallel 
to  the  acceleration  sensitivity  vector.  An  impor¬ 
tant  result  of  Eq.  1  is  that  the  frequency  shift, 
f(a)-fo.  is  zero  for  any  acceleration  in  the  plane 
normal  to  the  acceleration  sensitivity  vector. 

The  frequency  shift  described  in  Eq.  1  is  also 
induced  by  the  acceleration  due  to  gravity 
(even  without  motion).  This  is  commonly  re¬ 
ferred  to  as  "2g-tipover.'  During  '2g-tipover,' 
the  magnitude  of  the  acceleration  is  1g  in  the 


direction  towards  the  center  of  the  earth.  (The 
magnitude  of  acceleration  given  in  this  paper 
will  be  in  units  of  g,  i.e.,  the  magnitude  of  the 
earth's  gravitational  acceleration  at  sea  level, 
980  cm/sec^.)  When  a  resonator  is  rotated 
180“  about  a  horizontal  axis,  the  scalar  product 
of  the  acceleration  and  the  unit  vector  normal 
to  the  initial  lop”  of  the  resonator  goes  from 
•1g  to  -t-lg,  which  is  a  difference  of  2g. 

Fig.  1  shows  actual  data  of  the  fractional 
frequency  shift  of  a  resonator  (operating  in  an 
oscillator)  when  the  oscillator  is  rotated  about 
*hree  mutually  perpendicular  axes  in  the  earth's 
gravitational  field.  For  each  curve,  the  axis  of 
rotation  is  held  horizontal.  The  sinusoidal 
shape  of  each  curve  is  a  consequence  of  the 
scalar  product  being  proportional  to  the  cosine 
of  the  angle  between  the  acceleration  sensitivi¬ 
ty  vector  and  the  acceleration  due  to  gravity. 


Figure  1.  ”2-g  tipover,”  frequency  change 
versus  rotation  in  Earth's  gravitational  field  for 
three  mutually  perpendicular  axes. 


THE  EFFECT  OF  VIBRATION  ON  A  CRYSTAL 
RESONATOR 

In  nf\ost  applications,  the  magnitude  of  the 
acceleration  is  time-dependent.  The  magnitude 
of  the  acceleration  sensitivity  vector,  for  accel¬ 
eration  amplitudes  commonly  encountered,  is 
indeper\dent  of  acceleration  amplitude."  The 


time-dependent  trequency  shift  due  to  a  com¬ 
plex  vibration  can,  therefore,  be  determined 
from  the  sum  of  the  individual  sinusoidal  com¬ 
ponents,  i.e.,  the  system  is  linear  and  superpo¬ 
sition  holds.  Simple  sinusoidal  vibration  will  be 
discussed  first;  the  extension  will  then  be  made 
to  random  vibration. 

Simple  harmonic  motion  will  be  assumed,  with 
an  acceleration  given  by 

a  *  A  cos(2ji/^f), 

where  X  is  the  peak  acceleration  vector  in  units 
of  g,  is  the  frequency  of  vibration  in  Hertz, 
and  t  is  time  in  seconds.  The  variation  of  the 
frequency  with  time  can  be  determined  by 
combining  Eq.  2  with  Eq.  1 ,  resulting  in 

f(a)  »  4(1  A  )  cos(2ji4f)).  (3) 

The  behavior  of  the  device  can  be  described 
by  defining  a  rectangular  coordinate  system 
fixed  to  the  resonator.  The  acceleration  sensi¬ 
tivity  vector  and  the  acceleration  vector  can 
then  be  described  in  terms  of  the  three  unit 
vectors  defined  by  that  coordinate  system. 
Therefore  Eq.  3  can  be  transformed  into  a 
scalar  equation  containing  the  thiee  (i,  j,  and  k) 
components  of  X  and  f,  i.e., 

f(a  )-4*(1  *cos(2nf/). 

1  his  can  be  rewritten  as 

f(a)»4*Af  cos(2ji4f), 

where 

Af « (6) 

is  the  peak  frequency  shift  due  to  the  accelera¬ 
tion  K. 

There  a.m  three  quantities  with  units  of  frequen¬ 
cy  to  keep  in  mind:  4,  fy.  and  Af.  It  can  be 
seen  from  Eq.  5  that  the  output  frequency 
deviates  from  the  center  frequency,  4,  by  the 


amount  ±  Af,  at  a  rate  of  This  Is  shown 
schematically  In  Fig.  2. 
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Rgure  3.  (a)  Acceleration  level  versus  time, 
(b)  Resulting  c -dilator  output  showing  frequen¬ 
cy  modulation. 


Rgure  2.  Instantaneous  "carrier”  frequency  for 
several  instants  during  one  cycle  of  vibration. 

Each  plot  is  the  instantaneous  output  frequency 
of  a  crystal  osdilator  while  undergoing  a  vibra¬ 
tion  at  frequency  Rg.  3b  is  the  voltage 
versus  time  at  the  output  of  the  same  crystal 
osdilator  showing,  in  a  much  exaggerated  way. 
the  variation  in  the  frequency  with  acceleration 
amplitude.  Rg.  3a  is  the  acceleration  wave¬ 
form. 

FREQUENCY  DOMAIN 

It  is  very  useful  to  transform  the  effect  of  vibra- 
tiori  into  the  frequency  domain.  This  will  allow 
the  formulation  of  a  oonveniem  measurement 
scheme  and  allow  comparison  of  vibration 
effects  to  more  familar  forms  of  phase  noise. 
The  voltage  appearing  at  the  output  of  an 
osdilator  is  given  by 

V-WHeosWO).  (7) 


I 

♦(/).2icjf(/')df'.  (8) 

When  the  osdilator  frequency  is  varying  due  to 
simple  harmonic  acceleration  modulating  the 
resonant  frequency  of  the  resonator,  the  phase 
in  Eq.  7  becomes,  using  Eqs.  5  and  8. 

♦(f)  -2xf,f^(Af/fJsln(2jc/',f).  (9) 

When  Eq.  9  is  inserted  into  Eq.  7.  the  result  is 
V{f)-V,cos(2x4f^{Af/fJsin(2K/,f)).  (10) 

Eq.  1 0  is  the  expression  for  a  frequency-modu¬ 
lated  signal.  It  can  be  expanded  in  an  Infinite 
series  of  Bessel  functions'*  resulting  in'* 
V'(r)-V4IJ,(P)cos(2jt4f) 

♦  J,(P)cos(2x(4<)/) 

♦  J,0)cos(2x(/,-4)/)  (11) 

♦  j^(P)cos(2x(4-24)/) 

♦  j^(P)cos(2x(/,-24)f )-...) 


where  the  phase  ♦(t)  is  derived  from  the  fre¬ 
quency  by 


3 


where  p  -  Alf^  ■  (^  *  Ibe  modulation 

Index  (from  standard  FM  theory). 


» 


The  first  term  In  Eq.  1 1  Is  a  sine  wave  at  the 
carrier  frequency  with  an  amplitude,  relative  to 
Vg,  of  Jo(p).  The  other  terms  are  vibration- 
induced  sidebands  at  frequencies  ^  -i-  f^,  f,  -  f„, 
fg  +  2f^  fg  -  2f^  etc.  (The  "sidebands"  for 
sinusoidal  vibration  are  spectral  lines,  which 
are  also  referred  to  as  "brightlines".)  The  ratio 
of  the  power  in  the  nth  vibration-induced  sideb¬ 
and  to  the  power  In  the  earner,  denoted  by 
l„{U.  Is  given  by 

or,  more  commonly  expressed  in  decibels  as 
/.{O(dft^-20log{J,(P)/Jg(P)).  (13) 

SMALL  MODULATION  INDEX 

Several  approximations  can  be  made  if  the 
modulation  Index  Is  less  than  0.1.  This  is  the 
case  for  most  frequency  standards  In  the  HF 
band  experiencing  "normal"  accelerations  of  10 
g  or  less  at  acceleration  frequencies  above  a 
few  hertz.  The  approximations  are 

Jg(P)-1;  p<0.1 

J,(P)-p/2;  p<0.1  (14) 

JM’O;  p<0.1,n>2. 

Therefore,  after  combining  Eqs.  13  and  14: 

/,(fj- 20.log  ((f -A  )V(20) 

Li(0-o 

For  small  modulation  Index,  therefore,  most  of 
the  power  Is  In  the  carrier  and  a  small  amount 
is  in  the  first  sideband  pair  (upper  and  lower). 
The  amplitudes  of  the  higher  sidebands  are 
negligible. 

As  an  example,  consider  a  5  MHz  oscillator 
with  a  vibration  sensitivity  of  2x10'*/g.  If  the 
oscillator  experiences  an  acceleration  of  5  g 
along  the  direction  of  the  acceleration  sensitivi¬ 
ty  vector,  the  value  of  Af  Is  (5x10*  Hz)  *  (5  g)  * 
(2x10'*/g)  -  0.05  Hz.  The  modulation  index 


and  l,(fj  for  several  vibration  frequencies  are 
shown  In  Table  1. 


TABLE  1 


First  sideband  levels,  in  dBc,  for  a  5  MHz  oscil¬ 
lator  with  an  acceleration  sensitivity  of  2x10'‘/g 
experiencing  an  acceleration  of  5  g  peak. 


UHz) 

P 

l,{fJ{dBc) 

5 

0.01 

-46 

25 

0.002 

-60 

50 

0.001 

-66 

500 

0.0001 

-86 

It  can  be  seen  that  I,  falls  off  at  6  dB  per 
octave  (20  dB  per  decade.)  A  consequence  of 
linearity  at  low  vibration  levels  is  that  these 
sinusoidal  accelerations  are  Independent.  The 
amplitude  at  any  frequency  is  the  same  If  the 
brightlines  are  excited  separately  or  In  combi¬ 
nation.  The  first  three  of  these  vibration-in¬ 
duced  spectral  lines  and  the  carrier  are  shown 
In  a  typical  spectmm  analyzer  output  In  Fig.  4. 


Hgure  4.  Vibration-induced  spectral  fines  and 
carrier  resulting  frorri  'Inusoidal  acceleration  at 
5,  25,  and  50  Hz. 
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FREQUENCY  MULTIPLICATION 

Phase  noise  during  vibration  is  of  great  con¬ 
cern  In  systems  such  as  radar,  navigation,  and 
satellite  communications,  where  the  frequency 
of  the  crystal  oscillator  Is  multiplied  up  to  the 
microwave  region.  Upon  frequency  multiplica¬ 
tion  by  a  factor  of  N,  the  vibration  frequency  f„ 
is  unaffected  since  it  is  an  external  influence. 
The  peak  frequency  change  due  to  vibration, 
Af,  however,  becomes 

AN{f-A)/V4. 

The  modulation  index  ^  is  therefore  increased 
by  a  factor  of  N.  Expressed  in  decibels,  fre¬ 
quency  multiplication  by  a  factor  N  increases 
the  phase  noise  power  by  20  log  N. 

The  relationship  between  the  vibration-induced 
phase  noise  power  of  two  oscillators  with  the 
same  vibration  sensitivity  and  different  carrier 
frequencies  is 

(17) 

where  lA(fv)  is  the  sideband  level,  in  dBc,  of  the 
oscillator  at  frequency  f^,  and  Ia(U  is  the  side¬ 
band  level  of  the  oscillator  at  frequency  fg. 
Frequency  multiplication  to  a  higher  frequency 
is  indistinguishable  from  direct  frequency  gen¬ 
eration  at  the  higher  frequency.  For  example, 
when  a  2x10'*/g  sensitivity  5  MHz  oscillator's 
frequency  is  multiplied  by  a  factor  of  315  to 
generate  a  frequency  of  1575  MHz,  its  output 
will  have  identical  vibration-induced  sidebands 
to  a  1575  MHz  SAW  oscillator  with  a  sensitivity 
of  2x10‘*/g.  This  is  identical  to  the  "20  log  N" 
term  associated  with  the  increase  in  phase 
noise  due  to  frequency  multiplication.  Again, 
this  relationship  holds  only  if  p  <  0.1. 

LARGE  MODULATION  INDEX 

If  the  modulation  index  p  is  larger  than  0.1 ,  the 
approximations  made  in  Eq.  14  are  not  valid. 
This  often  occurs  in  UHF  and  higher  frequency 
systems.'  At  5  g  acceleration,  for  example,  the 


vibration-induced  sidebands  produced  by  the 
1575  MHz  oscillator  and  the  5  MHz  oscillator 
multiplied  by  315  are  shown  in  Table  2. 


TABLE  2 

First  sideband  levels,  in  dBc,  for  a  1575  MHz 
oscillator  with  an  acceleration  sensitivity  of 
2x10''/g  experiencing  a  peak  acceleration  of 
5g. 


UHz) 

P 

Ii(fv)(dBc) 

5.25 

3.15 

+2.3 

25 

0.63 

-9.6 

50 

0.315 

-16.0 

500 

0.0315 

-36.0 

It  can  be  seen  that  the  values  of  I,(fJ  in  Table 
2  are  much  greater  than  those  in  Table  1  for 
the  same  vibration  frequency.  This  is  a  conse¬ 
quence  of  the  ratio  of  carrier  frequency  to 
vibration  frequency  being  much  larger.  It  is 
possible,  as  can  be  seen  in  this  example,  for 
the  sidebands  to  be  larger  than  the  carrier. 
Also,  there  are  even  conditions  where  the 
carrier  disappears  and  the  value  of  l,(fy)  goes 
to  infinity,  e.g.,  when  p  equals  2.4,  that  is,  all  of 
the  power  is  in  the  sidebands  and  none  is  in 
the  carrier. 

ALLAN  VARIANCE 

The  effect  of  sinusoidal  phase  modulation  on 
the  two-sample  deviation  (also  called  the 
square-root  of  the  Allan  variance)  aJx)  of  a 
frequency  standard  has  been  shown''  to  be 

Oy  ■  (^  /Ji  4 t)  sin*  (n  f^x)  ( 1 8) 

where  O  is  the  peak  phase  deviation,  f^  is  the 
frequency  of  the  phase  modulation,  and  t  is  the 
measurement  time.  It  can  be  seen  from  Eq.  9 
that  the  mag  litude  of  the  maximum  phase 
deviation  for  a  single  vibration-induced  side¬ 
band  is  Af/fy.  Therefore,  the  vibration-induced 


I 


» 


I 


I 


I 


I 


I 


I 


» 


5 


I 


a,(T) 


» 


Oy(t)  Of  a  frequar^  standard  with  an  accelera¬ 
tion  sensitivity  of  f,  subjected  to  acceleration  K 
at  a  frequency  of  f,,  is 

Oy -I(f  A  )/(s/^x)J sin*  (jtCx) . 

This  effect  is  shown  in  Rg.  5.  The  frequency 
standard  is  assumed  to  have  an  Allan  variance 
of  (1x10'^’/x)*  when  not  being  accelerated. 


tion  can  be  described  by  its  power  spectral 
density  G(f).  The  power  spectral  density  of  fre¬ 
quency  fluctuationr  S«(f),  can  be  obtained  by 
multipi^ng  the  power  spiral  density  of  accel^ 
oration  by  the  square  of  the  acceleration  sensi¬ 
tivity  in  (Hz/g)*.  The  single-sideband  measure 
of  stability,  (f),  Is  related  to  S«(f).  if  the  modu¬ 
lation  index  Is  small,  by^* 

Sf(0-S,(0/(2/*).  (20) 


A  plausibiKty  argument  for  the  occurrence  of 
the  peaks  and  valleys  in  Rg.  5  is  that  the  aver¬ 
age  frequency,  as  given  in  Eq.  5.  is  fo  for  an 
averaging  time  equal  to  an  integer  multiple  of 
the  period  of  the  vibration.  Since  this  is  a  con¬ 
stant.  the  variance  is  zero.  The  average  fre¬ 
quency  departure  from  f^  is  a  maximum  when 
the  averaging  time  is  an  integer  multiple  of 
one-half  the  period  of  the  vibration.  Therefore, 
the  peaks  in  Rg.  5  occur  when  x  >  (2n 
1)/(2fJ  and  the  valleys  occur  when  x  >  (n 
1)/4.  where  n  >  0, 1, 2,  3,.... 


Therefore.  tf(f)  for  random  vibration  is  given  by 

a(0-(|riO' G(0/(2^*)- 

As  an  example,  consider  the  random  accelera¬ 
tion  spectrum  given  in  the  upper  right  of  Rg.  6. 
G(f).  in  units  of  g^Hz,  is  given  by 

G(/)-0.04  5</<220Hz 

G(/)-C.07*(f/300)*  220</<300Hz  ,22) 

G(/)-0.07  300</<1000Hz 

G{ /)- 0.07*  (//1 000)-*  1</<2kHz. 


Rgure  5.  o,(x)  versus  x  for  an  oscillator  with 
acceleration-sensitivity  of  1  x  10'*/g  during  1  g 
vibration  at  20  Hz. 


If  the  vribration  sensitivity  is  1x10'*/g  and  the 
oscillator  is  operating  at  10  MHz.  the  peak  fre¬ 
quency  deviation  is  0.01  Hz.  Therefore.  tf(f)  is 


2(0 


5<.<220Hr 
220<f< 300  Hr 

/(o.7KOJi)«j  300<  1 000  Hr 

(0jD7M0J1«1000)' 

i7^ 


]  1<f<2kHr 


which  is  shown  in  Rg.  6.  Note  that  outside  of 
the  vibration  frequency  range  defined  by  the 
given  acceleration  spectmm,  the  phase  noise 
is  identical  to  that  of  a  non-acceletated  devio). 


FW4DOM  VIBRATIONS 

In  most  situations  the  acceleration  experienced 
by  a  frequency  standard  is  not  simple  harmonic 
motion;  it  is  random,  i.e.,  the  vibratory  power  is 
randomly  otistributed  over  a  range  of  frequen¬ 
cies,  phases,  and  amplitudes  The  accelera- 


INTEGRATED  PHASE  NOISE  AND  PEAK: 
PHASE  EXCURSION 
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Spedafists  in  crystal  resonators  and  oscillators 
generally  characterize  phase  noise  by  S^(f)  or 
Sf(f).  Many  users  of  crystal  oscillators,  howev- 
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RguraG.  Random  vibration^nduced  phase 
noise  and  random  vibration  enveiope. 


Therefore. 

O, »  1.4  jf  lO"*  racSans.  (27) 


Whiie  the  osdilator  is  vibrating,  the  phase 
noise  is  given  by  Eq.  25  in  the  ba^  from  1  Hz 
to  5  Hz.  and  by  Eq.  23  in  the  band  from  5  Hz 
to  2000  Hz.  The  integrated  phase  noise  is 


*,• .  2  a .  Isa^df 

1  • 

900  1000 

♦  ji0.07)l!^fdf* 

tooe 

^  r  (Q.07K0.01«10C0)‘^ 

1000 

■  8  X  10'^  racSans^. 


(28) 


er,  characterize  phase  noise  in  terms  of  'phase 
jitter."  Phase  Jitter  is  the  phase  noise  inte¬ 
grated  over  the  system  bandwidth.  Similarly,  in 
phase-lock-ioops,  it  is  the  magnitude  of  the 
phase  excursions  that  determines  whether  or 
not  the  loop  will  break  lock  under  vibration. 

One  can  use  the  previous  example  to  investi¬ 
gate  the  effect  of  vibration  on  integrated  phase 
noise.  Integrated  phase  noise  is  defined,  for 
the  band  f,  to  f,.  as 

k 

•  j  SJif)df,  (24) 

i 

where  is  the  spectral  density  of  phase, 
equal  to  2  SL  in  the  frequency  band  of  1  Hz  to 
2  KHz,  the  phase  noise  of  the  nonvibrating 
oscillator  from  Rg.  6  is  given  by 

/si  kHz  (25) 

a»1x10***  /fclkHz 

and  the  integrated  phase  noise  in  the  same 
band  is 

-  2  X  20*’*  ra<Mans* .  (26) 


Therefore. 

■  9  X 10"*  ratSans.  (29) 

While  the  osdiiatnr  is  vibrating,  it  can  be  seen 
that  the  integrated  phase  noise  is  4000  times 
that  of  the  noise  when  it  is  not  vibrating  and 
the  rms  phase  deviation,  is  about  60  times 
larger  during  vibration. 

When  the  osdiiator  is  subjected  to  a  simple 
sinusoidal  vibration,  the  pe^  phase  excursion 
follows  from  Eq.  9,  i.e., 

0^-A///^.  (30) 

For  example,  if  our  10  MHz,  1x10'*/g  osdiiator 
is  subjected  to  a  10  Hz  sinusoidal  vibration  of 
amplitude  1  g,  the  peak  vibratbn-induced 
phase  excursion  is  1x10'’  radians.  If  this  osdi¬ 
iator  is  used  as  the  reference  osdiiator  in  a  10 
GHz  radar  system,  the  peak  phase  excursion 
at  10  GHz  will  be  1  radian.  Similarly,  when  an 
osditator's  frequency  is  multiplied,  the  square 
root  of  the  integrated  phase  noise  is  increased 
by  the  multipGcation  factor.  At  10  GHz,  for 
example,  ^  in  the  above  example  increases 
from  OxlO'*  radians  to  0.9  radiansi  Such  large 
phase  excursions  are  detrimental  to  the  perfor¬ 
mance  of  many  systems,  such  as  those  which 
employ  phase  lock  loops  or  phase  shift  keying. 
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MEASUREMENT 


The  sidebands  generated  by  vibration  can  be 
used  to  measure  the  acceieration  sensiti\^. 
Eq.  15  can  be  rearranged  to  get 


r,-(2/,/A,4)io««'”. 

(31) 

where  F,  is  the  component  of  the  acceieration 
sensitivity  vector  In  the  i  direction.  Three  mea¬ 
surements,  aiong  mutually  perpendicular  axes, 
are  required  to  characterize  f.  which  becomes 

(32) 

with  a  magnitude  of 

|f 

(33) 

One  scheme  for  measuring  T  is  shown  in  Fig. 
7.  The  local  oscillator  Is  used  to  mix  the  carrier 
frequency  down  to  the  range  of  the  spectrum 
analyzer.  If  the  local  oscillator  Is  not  modulat' 
ed,  the  relative  sideband  leveis  are  unchanged 


Hgure  7.  Typical  equipment  configuration  for 
measurement  of  accelerTtion-sensitivity  using 
sinusoidal  acceleration. 

by  mixing.  The  frequency  multiplier  is  used  to 
overcome  dynamic  range  imitations  of  the 
spectmm  arudyzer,  using  the  *20  log  N'  en- 
hancerrwnt  (fiscussed  previously.  The  mea¬ 


sured  sideband  ieveis  must  be  adjusted  for  the 
multiplication  factor  prior  to  insertion  into  Eq. 
31.  It  must  be  stressed  that  Eq.  17  Is  valid  only 
if  ^  <  0.1.  A  sampie  measurement  output  and 
caiculation  is  given  in  Rg.  8. 


Rgure  8.  Sample  acceleration-sensitivity  mea¬ 
surement  results  for  sinusoidal  acceleration. 


Other  measurement  schemes  Include  passive 
excitation’^  of  the  resonator,  the  use  of  random 
vibration,’*  and  the  2-g  tipover  tost.  Many 
oven-controlled  oscillators  are  not  suitable  for 
characterization  by  the  2-g  tipover  test,  hov.vev- 
er,  because  rotation  of  the  oscillator  results  in 
temperature  changes  Inside  the  oven  that  can 
mask  the  effects  due  to  the  acceleration  chang¬ 
es.  Vibration  tests  are  also  subject  to  pitfalls 
since  resonances  in  the  oscillator  or  shake 
table  assembly  can  produce  false  results.  It  is. 
therefore,  important  to  perform  the  test  at  more 
than  one  frequency. 


THEORY 

The  theoretical  understanding  of  acceleration 
sensitivity  does  not  yet  enable  one  to  predict 
the  acceleration  sensitivity  of  'rear  resonator 
designs.  Theoretical  activity  started  with  the 
study  of  in-plane  forces  on  simplified  resonator 
models.  The  first  studies  were  concerned  with 
the  dosely  related  force-frequency  effect.  This 
effect,  first  reported  by  Bottom’*  in  1947,  is  the 
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change  in  frequency  induced  by  a  pair  of  op* 
posed  forces  in  the  plane  of  a  resonator  plate 
applied  at  the  rim.  Those  forces  distort  the 
quartz  plate  and,  because  of  the  nonlinear 
elastic  behavior,  change  the  acoustic  velocity.'^ 
Since  the  frequency  of  a  resonator  is  a  function 
of  the  acoustic  velocity  and  the  dimensions  of 
the  quartz  plate,  the  forces  change  the  fre* 
quency. 

The  first  attempt  at  an  analytical  solution  to  the 
force-frequency  problem  was  made  by  Mingins, 
Barcus,  and  Perry  in  1962.“  They  assumed 
an  AT-cut  crystal  plate  of  infinite  lateral  dimen¬ 
sion  and  used  a  perturbation  technique  with 
linear  elastic  coefficients.  In  1963,  they  dis¬ 
cussed  nonlinear  theory  and  the  need  to  in¬ 
clude  the  third-order  elastic  coefficients.^' 
Those  coefficients  were  measured  in  1966  by 
Thurston,  McSkimin,  and  Andreatch.“  The  first 
calculation  to  use  the  nonlinear  theoiy  was  per¬ 
formed  by  Keyes  and  Blair  in  1967.“  In  1973, 
Lee.  Wang,  and  Markenscoff^*  published  the 
first  of  a  continuing  ssries  of  papers  (by  Lee 
and  his  students),  in  which  they  calculated  the 
force-frequency  coefficient  (as  defined  by 
Ratajski”)  as  a  function  of  azimuth,  using  the 
general  theory  oi  incremental  elastic  deforma¬ 
tions  superimposed  on  finite  initial  deforma¬ 
tions.  The  next  year,  lee  et  al  investigated  the 
frequency  change  due  to  cantilever  bending  of 
the  plate.“  In  1977,  a  variational  analysis  of 
the  force-frequency  effect,  including  the  effect 
of  material  anisotropy  on  initial  stress,  was 
given  by  EerNisse,  Ballato,  and  Lukaszek^^  for 
doubly  rotated  aits.  In  1978,  Janiaud,  Nissim, 
and  Gagnepain”  obtained  analytic  solutions  for 
the  biasing  stress  in  singly  and  doubly  rotated 
plates  subject  to  diametric  forces.  Those 
results  allowed  them  to  calculate  Ratajski's 
force-frequency  coefficient,  and  also  the  in¬ 
plane  acceleration  sensitivity.  Lee  and  Wu“ 
extended  the  work  of  reference  24  to  treat 
plates  of  any  cut.  In  that  paper,  the  solution  for 
an  arbitrary  number  of  ribbon  supports  was 
obtained. 


Under  acceleration,  body  forces  in  the  quartz 
plate  are  balanced  by  reaction  forces  from  the 
mounting  structure.  As  in  the  force-frequency 
effect,  distortion  of  the  crystal  lattice  causes  the 
resonant  frequency  to  change.  For  bulkwave 
resonators,  only  the  in-plane  case  has  been 
treated.  The  first  analysis  was  by  Valdois, 
Besson,  and  Gagnepain"  in  1974,  who  showed 
that  the  effect  is  linear.  Several  papers  by  Lee 
and  Wu“‘®’  considered  resonators  with  three- 
and  four-point  mounts.  Two-point  mounts  were 
considered  by  Janiaud,  Nissim,  and 
Gagnepain.“  For  SAW  resonators,  recent 
work  by  Shick  and  Tiersten^  and  by  Sinha  and 
Locke“  has  treated  acceleration  both  in  the 
plane  of  the  plate  and  normal  to  it. 

There  seems  to  be  a  significant  dependence  of 
the  acceleration  sensitivity  on  small  variations 
in  the  mounting.  Analyses  of  bulkwave  resona¬ 
tors,  to  date,  have  considered  only  point  sup¬ 
ports.  Real  resonators  have  supports  that  are 
distributed  over  a  finite  area.  In  addition,  initial 
stress  conditions  are  difficult  to  determine  and 
may  not  be  well-reproduced.  There  is  an 
ongoing  effort  by  Lee  and  Tang**  to  use  finite 
element  analysis  to  more  accurately  model  the 
mounting  structure.  It  is  hoped  that  this  will 
allow  computer  simulations  that  will  permit  the 
detemiination  of  the  optimum  design  parame¬ 
ters  for  minimizing  acceleration  sensitivity. 


EXPERIMENTAL  RESULTS 

Several  papers  have  reported  experimental 
results  on  the  force-frequency  effect""”  *  and 
the  effects  of  bending  moments.*  All  used 
point  mounts  in  special  fixtures,  and  results 
agree  fairly  well  with  the  theoretical  analyses. 
The  reported  experimental  observations  of  the 
actual  acceleration  sensitivity  of  real  resona¬ 
tors.  on  the  other  hand,  is  remarkable  in  that  it 
defies  simple  explanation. 

The  experimental  effort  on  acceleration  sensi¬ 
tivity  started  in  the  1960's  with  Smith,  Spencer, 


9 


and  Warner.’’®  Valdois,  Besson,  and 
Gagnepain”  made  measurements  on  resona¬ 
tors  and  oscillators  to  demonstrate  that  the 
resonator  was  the  acceleration-sensitive  ele¬ 
ment.  Results  from  2-g  tipover  experiments 
were  reported  for  SC-cut  resonators  by  Kus- 
ters,  Adams,  Yoshida,  and  Leach®^  in  1977.  In 
1979,  Warner,  Goldfrank,  Meirs,  and 
Rosenfeld”  reported  results  on  SC-  and  FC-cut 
devices.  In  1981 ,  the  first  of  a  series  of  papers 
intended  to  explore  a  wide  range  of  resonator 
fabrication  paiameters  was  published  by  Filler 
and  Vig.®®  In  that  same  year,  Nakazawa, 
Lukaszek,  and  Ballato®®  reported  constant- 
acceleration  results  that  suggested  an  acceler¬ 
ation  sensitivity  that  is  not  linear  with  accelera¬ 
tion  level.  A  significant  parameter  for  the 
reduction  of  acceleration  sensitivity  was  report¬ 
ed  by  Filler,  Kosinski,  and  Vig^’  in  1982.  As 
one  makes  plano-convex  and  biconvex  AT-cut 
resonators  flatter,  the  acceleration  sensitivity 
decreases.  It  must  be  noted,  however,  that 
other  effects,  possibly  due  to  insufficient  "mode 
shape  control"  (see  section  on  that  subject 
below),  caused  large  scatter  in  the  data.  A 
survey  of  a  large  number  of  resonator  parame¬ 
ters,  such  as  SC-cut  contour,  thickness,  drive 
level,  temperature,  and  angle  of  cut,  was  pub¬ 
lished  in  1983  by  Filler,  Kosins!<i,  and  Vig.*^ 
No  significant  coi  relations  were  found. 

A  paper,  by  Weglein®*  in  1984,  reported  ♦ho 
acceleration  sensitivity  of  VHF  AT-  and  SC-cut 
resonators.  All  of  the  resonators  in  that  study 
were  disassembled  after  the  acceleration  sensi¬ 
tivity  measurements  to  examine  fabrication 
details.  That  investigation,  which  seems  to 
summarize  the  eiicrts  to  date,  showed  "little 
dependence  on  any  recognizable  parameter," 
except  that  the  lowest  average  acceleration 
sensitivity  was  found  in  the  group  of  AT  units. 
The  large  spread  in  the  data,  typical  of  all 
experimental  observations  of  acceleration 
sensitivity,  seems  to  be  caused  by  a  combina¬ 
tion  of  sometimes  offsetting  subtle  effects, 
which  are  difficult  to  control  during  fabrication. 


REDUCTION  OF  THE  ACCELERATION  SEN¬ 
SITIVITY  OF  RESONATORS 

The  intrcduction  of  the  stress-compensated  cut 
of  quartz,  i.e.,  the  SC-cut,  was  accompanied  by 
the  widespread  expectation  that  the  SC-cut 
would  have  significantly  lower  acceleration 
sensitivity  than  other  cuts.  Unfortunately,  that 
expectation  has  not  been  realized.  The  lowest 
ao^eleration  sensitivity  achieved  with  AT-cut 
resonators®’'’®  equals  that  achieved  with  SC-cut 
resonators;  both  cuts  can  have  low  parts  in 
10’®/g  sensitivity. 

Efforts  to  reduce  the  sensitivity  of  individual 
resonators  to  the  effects  of  acceleration  have 
stressed  the  support  structure.  Lukaszek  and 
Ballato®®  proposed  a  plate  geometry  that  would 
assure  the  proper  support  configuration  to 
reduce  the  force-frequency  effect.  Besson, 
Gagnepain,  Janiaud,  and  Valdois®®  proposed  a 
support  structure  that  insured  symmetry  with 
the  median  plane  of  the  resonator  plate. 
Oebaisieux,  Aubry,  and  Groslambert®®  a:>d 
Aubry  and  Debaisieux®^  have  reported  results 
using  QAS  resonators,  a  variation  of  the  BVA,'* 
which  insures  symmetry  of  mount  as  wall  as 
accurate  mount  locations  for  reducing  the 
force-frequency  effect.  Their  results  show  a 
marked  reduction  in  the  scatter  of  the  mea¬ 
sured  acceleration  sensitivity. 

MODE  SHAPE  CONTROL 

Recent  theoretical  and  experimental  evidence 
by  EerNisse,  et  al.,®’®®  and  Tiersten  and  Zhou*’ 
indicates  that  the  major  variable  yet  to  be 
controlled  properly  is  the  shape  and  location  of 
the  mode  (active  region  of  vibration)  in  a  "real 
world"  resonator.  Theoretically,  it  can  be 
shown  that  the  frequency  shift  due  to  a  stress 
on  the  resonator,  such  as  that  caused  by 
acceleration,  is  given  by  the  integral  of  the 
third-order  elastic  constant  times  the  mode 
energy  density  times  the  strain  pattern  set  up 
by  the  acceleration.  Leaving  out  tensor  nota¬ 
tion,  for  simplicity. 
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A///-C,  j^GHX,Y7)  E{X,Y,Z)dV  (34) 

where  Af/f  is  fractional  frequency  shift,  G  is  the 
normalized  strain  distribution  associated  with 
the  mode  of  vibration  (G^  is  proportional  to  the 
strain  energy  density),  E  is  the  strain  distribu¬ 
tion  associated  with  the  stress  caused  by 
acceleration,  C,  is  the  third-order  elastic  con¬ 
stant  tensor,  V  is  the  volume  of  the  resonator, 
and  dV  -  dXdYdZ. 

EerNisse  showed  that  very  small  changes  in 
the  mode  shape  and  location  have  significant 
effects  on  the  integral,  i.e.,  on  the  acceleration 
sensitivity  of  resonators.  He  further  argued 
that,  because  of  material  imperfections  (such 
as  dislocations),  and  fabrication  imperfections 
(thickness  variations,  contouring  flaws,  elec¬ 
trode  mass  density  and  thickness  variations), 
the  demands  on  mode  shape  control  will  be 
impractical  to  accomplish  during  manufacturing 
without  some  ability  to  Irim"  the  mode  shape 
and  location  for  each  individual  resonator.  He 
was  able  to  produce  changes  in  the  accelera¬ 
tion  sensitivity  of  resonators  by  depositing 
platinum  on  selected  areas  of  resonators  (in 
order  to  change  the  mode  shape  and  location). 

ACCELERATION  COMPENSATION  OF  OS¬ 
CILLATORS 

The  lack  of  progress  in  reducing  the  accelera¬ 
tion  sensitivity  of  the  resonator  below  the  low 
parts  in  tO'^/g  level  has  spawned  several  tech¬ 
niques  for  compensation  of  the  effect.  There 
are  two  general  classes  of  compensation, 
passive  and  active.  The  first  compensation 
results  were  published  by  Gagriepain  and 
Walls  in  1977.”  They  used  the  passive  meth¬ 
od  of  mechanically  arranging  two  resonators 
such  that  the  components  of  the  acceleration 
sensitivities  normal  to  the  plates  were  anti¬ 
parallel.  The  resonators  were  electrically 
connected  in  series  in  a  single  oscillator. 
Przyjemski”  and  Emmons”  used  an  active 
technique.  They  sensed  the  acceleration 


magnitude  with  an  accelerometer  aligned  with 
the  direction  of  the  acceleration  sensitivity 
vector  of  the  resonator.  The  accelerometer 
signal  was  fed  into  a  tuning  circuit  in  the  oscil¬ 
lator  in  order  to  counter  the  acceleration-in¬ 
duced  frequency  changes.  A  limitation  of  this 
technique  is  the  requirement  on  the  linearity  of 
the  tuning  network  at  all  operating  points. 
Emmons  also  employed  the  dual-resonator 
technique  and  suggested  using  an  accelera¬ 
tion-sensitive  capacitor  in  the  tuning  network  of 
the  oscillator.  The  latter  technique  is  available 
as  an  option  on  a  commercial  cesium-beam 
frequency  standard. 

Ballato”  suggested  a  method  for  compensation 
in  all  directions  using  a  resonator  pair  made  of 
enantiomorphs.  He  argued  that  opposite-hand¬ 
edness  is  the  only  way  to  line  up  all  three  crys¬ 
tallographic  axes  to  to  antiparallel.  This  was 
extended  by  a  series  of  patents  by  Ballato  and 
Vig.”^  A  simplification  was  patented  by 
Filler,”  who  showed  that,  since  the  accelera¬ 
tion  sensitivity  has  vector  properties,  all  that  is 
required  is  that  the  vectors  be  aligned 
antiparallel,  independent  of  the  handedness  of 
the  quartz.  Vig  and  Walls”  extended  this  work 
by  suggesting  a  method  to  accommodate  reso¬ 
nators  with  different  acceleration  sensitivity 
magnitudes. 

Rosati  patented”  an  active  technique  of  com¬ 
pensation  which  was  further  developed  by 
Rosati  and  Filler.*'  This  method  makes  use  of 
the  polarization  effect  in  doubly  rotated  resona¬ 
tors,  i.e.,  that  the  resonant  frequency  of  a 
doubly  rotated  resonator  is  a  function  of  the 
voltage  applied  to  the  electrodes.”  If  one 
senses  the  acceleration  using  an  accelerome¬ 
ter  and  feeds  that  signal  with  appropriate 
amplification  and  phase  reversal  directly  to  the 
resonator  electrodes,  compensation  can  be 
achieved.  The  advantage  of  feeding  the  cor¬ 
rection  signal  to  the  resonator  electrodes  rattier 
than  a  varactor  is  that  the  polarization  effect 
has  superior  linearity.  When  the  correction 
signal  is  applied  to  a  varactor,  the  nonlinearity 


of  the  varactor  and  of  the  frequency-capaci¬ 
tance  function  causes  sidebands  at  harmonics 
of  the  vibration  frequency.  One  implementation 
of  this  technique  was  used  to  compensate  a 
rubidium  oscillator.”  A  refinement  was  pro¬ 
posed  by  Frerking.” 

CONCLUSION 

Vibration  effects  are  a  significant  problem  in 
modem  communication,  navigation,  and  radar 
systems.  Progress  has  been  made  in  under¬ 
standing  the  causes  of  acceleration  sensitivity 
but  achieving  less  than  a  few  parts  in  10'°  per 
g  has  been  elusive.  Mode  shape  control  may 
allow  significant  improvements  in  the  future.  A 
great  deal  of  effort  has  been  expended  on 
compensating  for  acceleration  sensitivity,  but 
more  work  is  needed  to  improve  the  level  and 
bandwidth  of  the  compensation. 

An  IEEE  committee  has  prepared  a  report 
which  reviews  The  Effects  of  Acceleration  on 
Precision  Frequency  Sources,*  and  provides 
guidelines  for  the  specification  and  testing  of 
oscillator  acceleration  sensitivities.”*® 
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